Draft version April 16, 2013 

Preprint typeset using L^T^X style cmulatcapj v. 5/2/11 



ON THE TIDAL DISSIPATION OF OBLIQUITY 

T.M. Rogers 

Department of Planetary Sciences, University of Arizona, Tucson, AZ, 85719 



o 

(N 

U 

Oh 
< 

in 



W 

Oh 
i 

o 

CO 



> 

OO 

^t 

o 
m 

> 

X 



D.N.C. Lin 

Astronomy and Astrophysics Department, University of California, Santa Cruz, CA 95064 and 
Kavli Institute for Astronomy and Astrophysics and School of Physics, Peking University, China 

Draft version April 16, 2013 

ABSTRACT 

We investigate tidal dissipation of obliquity in hot Jupiters. Assuming an initial random orientation 
of obliquity and parameters relevant to the observed population, the obliquity of hot Jupiters does 
not evolve to purely aligned systems. In fact, the obliquity evolves to either prograde, retrograde or 
90° orbits where the torque due to tidal perturbations vanishes. This distribution is incompatible 
with observations which show that hot jupiters around cool stars are generally aligned. This calls into 
question the viability of tidal dissipation as the mechanism for obliquity alignment of hot Jupiters 
around cool stars. 

Subject headings: internal gravity waves, angular momentum redistribution, extra-solar planets, hot 
jupiters 



1. INTRODUCTION 

Hot Jupiters (with masses comparable to that of 
Jupiter and orbital periods less than a week or so) are 
found around 1-2 % of solar type stars. A widely adopted 
scenario for their origin is that these planets formed at 
several AU from their host stars and underwent inward 
migration due to tidal interaction with their natal disks 
(Lin et al. 1996). Another class of dynamical models 
assumes the hot Jupiters were relocated to close proxim- 
ity to their host stars through close encounters between 
planets, s ecular chaos, or Kozai resonance with compan- 
ion stars dRasio fc F ord 199 6|lWufc Mur ray 2003; Fab 
rycky & Tremaine 2007; Wu & Lithwick 2011; Naoz et al. 
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New observations are continually testing and con- 
straining these theories. To date the obliquity between 
planets and their host stars have been measured us- 
ing the Rossiter-McLaughlin effect in more than 50 sys- 
tems. Several hot Jupiters have been observed to have 
high obliquity, O (including retrograde orbits in which 
O > 7r/2), while many show alignment. In general, 
hot Jupiters around cool stars with effective temperature 
T* < 6, 250K, tend to be aligned, while those around hot 



stars, T* > 6, 250K, appear to be misaligned ( Winn et al. 
2010 1 . In order to accoun t for this dichotomy, Winnet al. 
12010) and subsequently, Albrecht et al. (2012) have sug- 



gested that 1) all hot Jupiters were relocated to close 
proximity to their host stars by one of the dynamical 
processes described above, resulting in a random distri- 
bution of obliquity, 2) the obliquities of those planets 
around cool stars are damped by efficient tidal dissipa- 
tion in the convective envelopes of these cool stars and 
3) the obliquities of those planets around hot stars, with 
radiative envelopes, reflects that of the initial random 
obliquity distribution because tidal dissipation in these 
systems is inefficient. 



Quantitatively, Albrecht et al. (2012) evaluated the 
magnitude of the tidal dissipation time scale in stars 
with convective envelopes, tce, and those with radia- 
tive envelopes, tra using models of equilibrium tides for 
convective solar-type stars an d dynamical tides for radia- 
tive, intermediate-mass stars ( Zahn|[l977 ). They showed 
a correlation between the magnitude of misalignment, O, 
and the dissipation timescales. 

Despite this suggestive evidence, there are large un- 
certainties i n the equilibrium-tide model that |Albrecht| 



et al. (2012) have adopted, mostly due to a controversial 



prescription for turbulent dissipatio n of the tides ( Gol- 
dreich &: Nichol son 1989; T erquem et al.|1998|[Ogilvie fe 



Lin||2007[|Barker fc Ogilvie||2010p . Using Zahn's model. 



Albrecht et al.| ( J2012 ) calculate that the magnitudes of 
Tce/t* and tra/t* vary by several orders of magnitude 
(see their Figures 24 & 25). The relatively small value 
of tce/t* obtained for many of the aligned hot Jupiters 
around solar type stars pose a challenge to their reten- 
tion against substantial orbital decay, unless the time 
scale for obliquity damping is substantially smaller than 
that for orbital decay. 
In a thorough theoretical analysis, 



Lai (|2012b 

a 



showed 
that the timescale for obliquity damping could be sub- 
stantially different from that of orbital decay. He iden- 
tified a component of the tidal torque which affects the 
obliquity alignment but not the orbital decay. He showed 
that, for the forcing frequency of this particular compo- 
nent of the tidal perturbation, inertial waves may be ex- 
cited to provide a dynamical tidal response which could 
lead to much more efficient energy dissipation and hence , 
shorter dissipation timescales. Therefore, the Lai (20121 
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theory appears to provide a potential mechanism to ac- 
count for hot Jupiters whose obliquity has been damped 
in absence of substantial orbital decay. 

In §2, we briefly recapitulate Lai's theory and in §3 
show that the orbital decay paradox can only be resolved 
if the timescale for obliquity evolution is much smaller 
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than that for semi-major axis decay. In the limit of neg- 
ligible orbital decay, we show in §3, that tidal dissipation 
does not lead to obliquity align ment . Bot h of these ef- 
fects have been clearly stated by Lai (2012 1. Our contri- 
bution is to provide the results of numerical integration 
for various limiting cases and compare them with the 
observational data. We discuss the implications of our 
results in §4. 

2. SEMI MAJOR AXIS DECAY AND OBLIQUITY 
ALIGNMENT 

We first briefly recapitulate the theory of equilibrium 
and dynamical tides raised by hot Jupiters on their host 
stars. 

2.1. Components of planets' tidal perturbation on their 
spinning host stars 

The tidal potential U(r,t) imposed by a planet with 
mass M p , with a circular orbit and an orbital angular 
frequency ft p on a star with a mass M* and a uniform 
spin with angular frequency ft* can be approximated to 
the lowest order in terms of spherical harmonics in a 
frame centered on the star with the z-axis parallel to the 
stellar spin such that 

U(r,t) =-J2 U mm ,(M p ,e)r 2 Y 2m (9,(f>)exp{-im'n p t) 

mm' 

(i) 

where the obliquity O is the angle between the stellar 
spin S and planet's orbital angular momentum vector 
L, whereas the magnitudes are given by S = /ft* and 
L = M p a 2 il p , I = kM*R 2 , where i?* is the stellar ra- 
dius, A;(~ 0.1) and I is the moment of inertia. In a frame 
co-rotating with the stellar spin, the forcing frequency is 
w mm ' = m'VLp — toSI* with seven components cont ribut- 



e. 



i ng to obliqu ity evolution ( Barker & Ogilvie 2009 



Lai ( |2012[ ) pointed out that it is possible lor some com- 
ponents ol the planets' tidal perturbing potential to have 
sufficiently small io mrn i < 217* to allow the corilois effect 
to provide the nec essary restoring fo rce for the excitation 
of inertial waves (Greenspan 1968). In cool stars these 
waves lead to dynamical tides ( Ugilvie fc Lin||2007 ) . For 
some forcing frequencies, the inertial waves ma y converge 
onto attra c tors and d issipate more efficiently ( Ogilvie & 
Lin||2004[ ) . |Lai| ( |2012[ ) showed that in the limit oi small 
ft*(< ft p ) , the only component of the perturbing poten- 
tial with sufficiently small uj mm > is that associated with 
(m,m!) = (1,0). For this component, tidal dissipation 
can lead to a shorter dissipation timescale and obliquity 
evolution without significant orbital decay. For more 
rapidly spinning stars (with ft* > ftp/2), it is possible 
to excite other components of the tidal response which 
would lead to both obliquity evolution and orbital decay. 

2.2. Evolutionary equations 

We first consider the possibility that the dissipation 
timescale, t mm > is indentical for all (to, rnf) components. 
In this case, dissipation of the equilibrium tide leads to 
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where r e = (Q Q /3fc 2 )(M*/Mj,)(a/i?*) 5 (P/27r) is the 
characteristic orbital evolution time scale, k 2 is the Love 
number, and Q a is the highly uncertain quality factor. 
In addition to turbulent dissipation of equilibrium tides, 
Q a may also include contributions from the dissipation 
of internal gravity a nd inertial waves in solar type stars 
(Ogilvie & Lin 2007). For rapidly spinning stars, the dy- 
namical response associated with the inertial waves may 
lead to the evolution of a and O on similar time scales. 
However, for slowly spinning stars (with ft* < ft p /2), 
only the forcing frequency associated with the (to, m ) = 
(1,0) component can excite the inertial modes, leading 
to a dynamical response, and relatively short timescales 
r io(<< T e)- The dynamical tide's contribution associ- 
ated with the (1,0) component of the torque leads to 
rates of change of a, ft*, and 9 such that 
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such that T W /r e - (Qw/Qa)(k 2 /k w )(S/L). Lai (2012) 
pointed out that since the (m,m') = (1,0) component 
does not contribute to a, obliquity alignment can oc- 
cur, in principle, prior to any significant orbital decay if 

no << t c . 



3. COMPUTATIONAL RESULTS 

3.1. Evolution of Semi-Major axis and Obliquity due to 
Equilibrium Tide 

Here we show that the obliquity alignment due to the 
equilibrium tide is accompanined by substantial orbital 
evolution. We first neglect any extra contribution from 
the (to, to') = (L0) component by numerically integrat- 
ing Equations pH), with a fourth order Runge-Kutta 
scheme, for an mital value of 9 = 45°. If we fix 
fc(i?*/a ) 2 (M*/Af p ) = 1, we have a set of solutions which 
depend only on the initial value of S/L. The solutions to 

(2) these integrations are shown in Figure 1. In that figure, 
black lines represent the evolution of a(t)/a , red lines 
represent the evolution of 9 and the blue lines represent 

(3) the evolution of the equilibrium timescale (which varies 
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as a~ 13 / 2 ) all as a function of time, in units of the (initial) 
equilibrium timescale. 

Concentrating first on the slow rotators (S/L <1) we 
see two important features. First, in order to have signif- 
icant obliquity damping the semi-major axis is reduced 
substantially. Given that we set our initial ratio {R„/a ) 
to 0.1 any reduction of a below this value implies the 
planet has fallen into the star. We see that this oc- 
curs at approximately 2.5 T eq o, a time when there is still 
non-negligible obliquity. Another way to look at this is 
to consider the time when the obliquity is damped to 
half its original value of tt/4, at this time the equilib- 
rium timescale has been reduced by nearly five orders of 
magnitude. One might hypothesize that we have caught 
the aligned hot Jupiters in their last gasp on their death 
march to infall in which their obliquity has been reduced 
to zero but their orbits have not been completely ex- 
hausted. However, the rapid decay of the equilibrium 
timescale coincident with this evolution is not consistent 
with the large population of aligned, yet surviving, sys- 
tems that have been observed. In summary, even modest 
amounts of obliquity damping are concurrent with signifi- 
cant orbital decay and a rapidly decreasing timescale over 
which that orbital decay will occur. Similarly, for fast 
rotators modest obliquity damping is accompanied by 
significant orbital expansion. Therefore, in order to ac- 
comodate any substantial obliquity evolution the planet 
would have to have started its orbit virtually inside the 
star, an impossible scenario. 

Therefore, for both slow and fast rotators it is impos- 
sible to get any substantial obliquity evolution without 
either 1) the planet falling into the star (slow rotators) 
or 2) the planet starting its orbital evolution inside the 
star (fast rotators). This long recognized problem leads 
to the conclusion that significant obliquity damping can 
only occur if the dynamical tide is considered and the 
timescale for obliquity alignment, no is substantially less 
than that for orbital decay/expansion, r eq . 

3.2. Population synthesis of obliquity evolution due to 
dynamical tides 



The model of Albrecht et al. (20121 assumes that hot 
jupiters arrived at their current, close positions with a 
random distribution of obliquities as a result of one of 
the dynamical processes listed in §1. To mimic this sce- 
nario we start from a random distribution of obliquities 
and integrate Equation (7), assuming a constant semi- 
major axis, a (this assumption is justified in the limit 
of no << T e ). Since the physics of tidal dissipation is 
highly uncertain we run a host of models varying three 
free parameters: S/L, rio/r eq and fl*/fl p . As stated 
above, in order for obliquity evolution to occur in the 
absence of orbital decay no << T e, so we consider only 
models for which no/ r e < 1 (although see comment be- 
low). We consider £l ast /£l p between 0.1 and 10. Finally, 
we consider values for S/L which vary between 0.1 and 
2. Note that these are all initial values. 

Figure 2 shows the results of one set of our integra- 
tions (see Figure caption for more details). The tidal 
potential generated by the planet on the star results in a 
torque which depends explicitly on the angle of obliquity 
and goes to zero if the obliquity is 7r, 7r/2 or 0. There 
are two components of this torque, one in the direction 
of the spin axis of the star, the other along the orbital 



axis of the planet. Which of these components dominates 
the evolution determines which eventual state of tt,tt/2 
or the planet tends to. Hot Jupiters with o < ir/2 
evolve toward alignment regardless the value of S/L after 
t ~ 2 — 10t 10 . On the other hand, most hot Jupiters with 
©o > 7r /2 either evolve towards it (in the limit of small 
S/L) or tt/2 (if S > L) after t ~ 2 - 10no- Therefore, 
if all hot Jupiters started with a random obliquity an- 
gle, tidal evolution would lead to a nearly equal division 
between those with prograde (O < 7r/2) and retrograde 
(O > 7r/2) or 90° orbits. The only circumstances under 
which all hot Jupiters evolve to aligned systems are when 
Tio > 0.5r e . However, as discussed in Section 3.1 such 
timescales would also lead to significant orbital decay and 
loss of planets to their host stars, which is inconsistent 
with the observations. 

3.3. Comparison with observation 

We obtain information from the website exoplanet.org 
and plot the distribution of obliquity as a function of 
S/L, fl*/Qp and stellar temperature in Figure 3. Plan- 
ets around solar type stars are represented by colored 
circles, while black circles represent planets around hot 
stars. A review of these orbital properties indicates that 
very few cool stars have S/L larger than 1. For these pa- 
rameters the most likely outcome of tidal dissipation of 
obliquity, is that initially random distributions of obliq- 
uity evolve to aligned, anti-aligned or 90° orbits. For 
S/L ~ 0.1 — 0.5 the percentage of mis-aligned systems 
is ~ 25 — 50%. The observations, on the other hand, 
show an overwhelming majority of hot Jupiters around 
cool stars have aligned orbits. The obliquities plotted 
in Figure 3 are of the projected obliquity and the error 
bars denote the likely range of true obliquity. Despite 
projection effects, the observations are still inconsistent 
with the theoretical prediction. 

The two limits considered, including the equilibrium 
tide or only the (1,0) dynamical tidal, bracket the so- 
lutions. We have avoided considering both components 
simultaneously as that would require assumptions about 
the efficiency of each mechanism. However, consider- 
ing the results of 3.2, that end states of obliquity evo- 
lution include a signficant fraction of retrograde orbits, 
and Equation (2), that indicates retrograde and prograde 
orbits may evolve on different timescales, one might ar- 
gue that retrograde systems might be preferentially lost. 
By inspection one can see that this depends on the ratio 
of fl^/tlp. We integrated Equations (2)-(4) varying the 
initial value of £1* /£l pi the results of which are shown in 
Figure 4. For f2*/f2 p < 1, where the majority of systems 
live, prograde and retrograde orbits evolve on similar, or 
at least indiscernible, timescales. For f^/fLj ~ 1 — 2, 
retrograde systems could be lost, while prograde sys- 
tems remain (with high obliquity), however this repre- 
sents a small fraction of the observed systems. Finally, 
for f2*/rj p >> 1 retrograde orbits are lost while prograde 
orbits migrate outward rapidly, inconsistent with the ob- 
served population of close-in, aligned systems. There- 
fore, we conclude that obliquity alignment due to tidal 
dissipation is unable to explain current observations. 

4. DISCUSSION 

The purpose of this paper is not to rule out the forma- 
tion of some hot Jupiter's through dynamical processes. 



Rogers & Lin 



The broad distribution of planetary eccentricity may be 
due to dynamical relaxation processes which can, in prin- 
ciple, lead to the scattering of a few gas giant planets to 
the proximity of their host stars. Kozai effect may also 
deliver gas giants, at least in some well known systems 
such as HD80606. Some of these encounters may be suf- 
ficiently close to induce strong tidal interaction between 
the planets and their host stars. It is entirely possible 
that dynamical processes may have led to the formation 
of some hot Jupiter with diverse obliquities around both 
solar type and hotter main sequence stars. 

However, the dynamical model alone cannot account 
for the dichotomy between the obliquity distribution of 
hot Jupiters around solar type and intermediate-mass 
stars because this mec hanism does not depe nd on the 
mass of the host stars. Albrecht et al. ( 2012| ) attribute 
this difference to the efficient dissipation of hot Jupiters' 
tidal perturbation on their solar type host stars whereas 
that process is inefficient in intermediate mass stars. The 
results presented here pose a challenge to this tidal align- 
ment model. 

We have shown that if one considers only the equilib- 
rium tide dissipation leads to such severe orbital decay 
that any systems which have their obliquity aligned fall 
into their host star, or, for fast rotators, will have their 
orbits expand so drastically that they would have to have 
started their orbit inside the star. At the other extreme, 
if we inste ad consider the dynamical tide proposed by 
Lai (2012) tidal dissipation of obliquity produces both 
prograde and retrograde systems for a population of hot 
Jupiters with random initial 0. We confirmed that al- 
though hot Jupiters with initially prograde obliquities 
(with 9 < 7r/2) would become aligned, those with ini- 
tially retrograde obliquities (with 8 > ir/2) would at- 



tain either anti-aligned or orthogonal obliquity. Such 
an obliquity distribution is inconsistent with that ob- 
served for hot Jupiters around solar type stars. We found 
that the possibility that these systems evolve by the dy- 
namical tide but that retrograde orbits are preferentially 
damped can only explain systems in which Q*/Q p ~ 1 
and there are few of these systems. Therefore, this also, 
can not explain the current observations. 

In a series of papers (Rogers et al 2012, 2013), we pro- 
posed an alternative model for the observed dichotomy 
between the spin-orbit alignment for hot Jupiters around 
solar-type and intermediate mass stars. We suggest that 
most hot Jupiters migrated to the proximity of their host 
stars through type II migration (Lin ei al. 1996). These 
planets retained the angular momentum vector associ- 
ated with the disk. But the spins of their hot host stars 
may be modulated due to the excitation, propagation, 
and dissipation of internal gravity waves, a process which 
is only efficient in hot stars. To date our gravity wave 
model is the only scenario which can provide a natural 
explanation of the observed difference between the dis- 
tribution of hot Jupiters around solar type stars and that 
around intermediate-mass stars. Furthermore, with this 
theory, there is no need to introduce multiple scenarios 
to account for the migration of hot Jupit ers versus that 
of mul tiple-planet systems as suggested by |Albrecht et al. 
((20131. 
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Fig. 1. — Evolution of semi-major axis, a (black lines), and obliq- 
uity (red lines) as a function of time for the equilibrium tide. 
Left axis represents semi-major axis and obliquity, right axis rep- 
resents equilibrium timescale. Left hand panel show slow rotators 
(S/L<1), while the right hand panel shows fast rotators (S/L>1). 
was initially set to 45° in these integrations. Different linetypes 
represent different values of S/L. 
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Fig. 2. — Obliquity evolution of a population of objects with initially random obliquity and with T-io/r eq = 0.001. Left panels show the 
original random distribution of obliquities for the 50 objects (open circles), while filled circles show the distribution of obliquity after 30r eq . 
Middle column shows the time evolution of the obliquity, while the right panel shows the fraction of objects as a function of time with 
prograde orbits (solid line), retrograde orbits (dotted line) and 90° orbits (dashed line). In all scenarios the objects evolve to prograde, 
retrograde or 90° orbits. Note that since we are keeping a fixed, S/L does not evolve. 
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Fig. 3. — Distribution of projected obliquity as a function of temperature, S/L and €l*/Cl p . Cool stars are represented in color, while hot 
stars are represented in black. 
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Fig. 4. — Evolution of obliquity and semi-major axis for various initial values of Q„/Clp, after 2T eq o. Red, blue and green lines represent 
initially retrograde, prograde and 90° orbits, respectively 



